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The  usual  mathematical  formulation  of  availability 
assumes  an  exponential  distribution  for  failure  and  repair 
tines.  While  this  is  sometimes  correct  for  reliability,  if 
is  not  likely  to  be  for  maintainability.  This  study  was 
conducted  to  verify  that  the  lognoraal  and  gaaaa  distribu¬ 
tions  are  suitable  descriptors  for  corrective  maintenance 
repair  times,  and  to  estiaate  the  difference  caused  in 
assuming  an  exponential  distribution  for  availability  and 
maintainability  calculations  when  in  fact  the  distribution 
is  lognoraal.  Port y- six  sets  of  data  of  electronic  and 
mechanical  systems  and  equipments  were  analyzed  using  the 
methods  of  probability  plotting  and  statistical  testing  for 
distributional  assumptions. 
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I . INTRODOCTIOH 


A.  OVERVIEW 

The  effectiveness  of  a  system  depends  not  only  on  its 
ability  to  aeet  its  specified  perfonance  requirements,  but 
also  on  its  ability  to  perform  when  needed,  for  the  duration 
of  its  assigned  missions,  and  for  its  operational  lifetime. 
The  technical  disciplines  concerned  with  these  time-related 
system  characteristics  are  reliability,  maintainability  and 
logistic  support. 

Commonly  used  methods  for  prediction  and  demonstration 
of  corrective  maintenance  repair  tines  assume  the  validity 
of  the  lognormal  distribution. 

In  order  to  help  in  focusing  attention  on  these  matters, 
a  statistical  analysis  on  data  sets  of  demonstrated  and 
field  repair  tines  has  been  conducted  as  part  of  a  study  on 
the  application  of  the  lognormal  and  gamma  distributions  to 
corrective  maintenance  data.  Ihe  results  of  the  analysis 
are  given  in  this  thesi3.  A  previous  study  of  electronic 
systems  and  equipments  concluded  that  the  lognormal  distri¬ 
bution  was  a  suitable  descriptor  for  repair  times  [Ref.  1]. 

B.  PURPOSE  AMD  APPROACH 

i •  Sbissiiisa 

The  objectives  of  this  study  have  been; 


(a)  To  test  the  exponential,  lognormal,  and  gamma 
distributions  as  descriptors  of  repair  time  data  on  mechani¬ 
cal  and  electronic  systems  and  equipments, 

(b)  to  verify  that  the  lognormal  or  gamma  distri  - 
but  ion  is  a  suitable  descriptor  for  repair  times, 

(c)  to  estimate  the  percentage  error  caused  in 
assuming  exponential ty  for  availability  and  maintainability 
calculations  when  in  fact  the  distribution  is  lognormal. 

2.  Systems  ajd  Data  Analyzed 

Forthy-six  sets  of  repair  time  data  for  electronic 
and  mechanical  systems  /  equipments  were  analyzed  (Table  1). 
Data  for  the  electronic  systems  (Sets  8-20)  came  from  formal 
maintanability  demonstration  test  reports  furnished  to  us  by 
the  Rone  Air  Development  Center,  O.S.  Air  Force.  Mechanical 
equipments  included  field  repair  time  pumps  used  in  French 
nuclear  electric  power  stations  (Sets  1-7),  and  additional 
26  sets  of  field  repair  data  from  British  Aircraft,  provided 
by  the  Onited  Kingdom  Ministry  of  Defence,  representing 
elapsed  tines  (Sets  21-33),  and  man  hours  (Sets  34-46).  So 
detailed  reports  are  available  for  the  mechanical  i-ems. 

3*  AaallS 1§  4£££°3£h 

Three  approaches  were  used  for  testing  data. 


Table  i 


SYSTEMS  /  EQUIPMENTS  ANALYZED 


Set  NO 


2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 
15a 

16 

17 

18 

19 

20 

21.34 

22.35 
23  \  36 

24 . 37 

25 . 38 
26 '39 
27;  40 

28.41 

29.42 
30;43 
3i;44 

32.45 

33 . 46 


SYSTEM  NAME 

Fr encE"Fe3s5nE5II"PuIps“ESwn"Time 

Prench  Pessenheim  Pomps  Repair  Time 

Condenser  Extraction  Pump 

Heater  Dram  Poap 

Main  Feedwater  Puip 

Reactor  Coolant  Poap 

Condenser  Vacuum  Puap 

AN/GRT-2  0 

AN/ARC-162  (?)  Radio  Set  Inter  mediate 
AN/ARC-163  (?)  UHP  Modern  Radio  Set 
AN/ALQ-1 25  TEREC  .  Direction  Finder  Set 
CGT  MAX.  HPU  and  R SS  Terminal  Organize 
AN/TPN-28  Radar  Beacon  Transponder 
AN/HST-I1  Electronic  Warfare  Training  Set 
AN/UYK-1  4C 
Revised  AN/UYK-1 4C 
Integrated  Microwave  G/G  Subsystem 
CCG-NBOr  Compass  Ears 
AN/GXS-2  7 
AN/ARN-1  01  (?) 

AN/FCC-9 8  Multiplexer 
Puselaga  Structure 
Plying  Controls 
Landing  Gear 
Hydra  ulics 

Ice  and  Ram  Protection 
Power  Plant  , 

Engine  starting 

Main  Rotor  Head  and  Blades 

Main  Rotor  Drive  ,  , , 

Electrical  Power  Supply  and  Distribution 

Co  maunications 

Surveillance  and  Search 

Automatic  Plight  Control 


Note:  In  sets  21-45,  the  two  numbers  represent  two 
part  data.  The  first  part  is  ELAPSED  TIME  and  second  part 
is  HAN  HOURS. 


(a)  Rough  check  using  a  technique  called  Boxplot. 
This  gives  a  rough  idea  that  the  data  cones 
from  symmetric  (normal)  or  skewed  distributions 
lognormal,  exponential  or  gamma). 

(b)  Plotting  the  data  using  appropriate  probability 
paper. 

(c)  Statistical  tests  of  significance  for  assumed 
distributions. 

The  following  procedure  was  used  for  analyzing  the 

data, 

i.  Sketch  boxplot  and  histogram  of  data  and 
determine  if  data  appears  to  come  from  a  symmetric  or  skewed 
distribution.  Obtain  an  idea  about  center  and  dispersion  of 
data. 

ii.  Plot  the  data  on  appropriate  probability 
paper  for  the  exponential  and  lognormal  distributions.  In 
some  cases,  plots  were  made  for  the  gamma  distribution. 

iii.  A  Chi-Square  Goodness-Of-Fit  test  was  per  - 
formed  for  lognormal  and  exponential  assumptions,  [Ref.  2] 

iv.  A  H  test,  due  to  Shapiro  and  »ilk,  [Ref.  3] 
was  used  to  test  the  lognormal  assumption  for  sample  sizes 
less  than  51,  (due  to  availability  of  the  tables) , 
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II.STAUgTICii  CONSIDERATIONS 
A.  ROUGH  CHECK  HETBODS 

1.  Boxplot  C  •  5  a  a  4  6] 

Purpose:  Many  batches  of  data  pile  up  in  the  middle. 
To  analyze  the  behavior  of  a  batch,  we  need  a  picture  of 
where  the  middle  (median)  lies  and  how  the  tails  relate  to 
it.  The  middle  is  generally  better  defined  than  the  rails 
so  we  want  to  see  more  at  the  tails.  Some  values  called 
outliers  are  so  low  or  so  high  that  they  seem  to  stand  apart 
from  the  rest  of  the  batch. 

One  technique  for  displaying  such  batches  of  data 
which  is  often  wore  convenient  than  a  histogram  is  the  box- 
plot.  But  the  histogram  is  still  good  for  showing  the  shape 
of  the  distribution.  Boxplot  will  not  show  multiple  modes. 
There  may  be  several  "modes",  and  we  want  to  see  them. 

A  boxplot  is  obtained  by  first  calculating  the  lower 
and  upper  quartiles  and  the  median  of  the  batch  (sample)  of 
numbers  and  then  plotting  tha  numbers  cn  a  horizontal  line. 

The  lower  quartile  is  the  value  that  divides  the 
batch  into  two  parts,  with  1/4  of  the  numbers  below  this 
value  and  3/4  above  it.  The  upper  qaartile  is  the  value 
with  3/4  of  the  observations  below  and  1/4  observations 
above  it. 
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The  next  step  is  to  draw  a  narrow  rectangular  box 


with  ends  corresponding  to  the  lower  and  upper  guartiles, 
and  to  display  the  median  point  by  a  plus  sign.  The  length 
of  this  box  is  called  the  interquartile  range  or  H-Spread 
(Figure  1)  . 

To  the  above  are  added  lines  marking  1/2  H-Spread  on 
each  side  of  the  box.  Data  values  outside  the  lines  (out¬ 
liers)  are  marked  with  stars  and  rectangular  circles. 

*  **  - 1  +  i -  *  *  0000 

Interquartile  range 
\ . 5  interquartile  *  range 

i  Interquartile  range  1 

-Figure  1  -  Boxplot- 

First,  the  data  are  ordered  such  that, 

X(1)  <  X  (2) .  <  X(n)  where, 

X ( 1 )  :  smallest  (minimum)  value 

X(n)  :  largest  (maximum)  value 

medians!!  ,  center  of  batch 

if  n  is  an  odd  number,  B=X((n  +  1)/2) 

if  n  is  an  even  number  ,  B»  (X  ( (N/2) +  1) +X  ( (n/2) -1) ) /2 
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guartiles=H  ,  center  of  half  of  batch 
d  (H)  *  (|  d  (M)  | ♦  1)  /2  where  1  A|  aeans  integer  part  of  & 
interguartile  distance  *  upper  quartile  -lower  guartile 

Here  are  good  boxplot  examples  for  the  exponential 
and  lognonal  model  (data  generated  randomly). 


BOZPLOT  EXPONENTIAL 


-I  ♦ 


-Figure  2  -  Boxplot  Exponential- 
BOZPLOT  NORMAL  (log  is  taken  then  boxplot) 


-Figure  3  -  Boxplot  Normal- 


2.  siatasia.  CM£-  5  6] 

Boxplot  is  good  for  comparison  but  a  histogram  is 
still  good  for  showing  shape.  Sometimes  histogram  is  more 


helpful  than  the  boxplot ,  especially  for  detection  of  more 
bumps  in  the  data.  Histogram  prints  data  in  classes.  Exam¬ 
ples  are  given  below  in  Figure  4. 

HISTOGRAM  EXPONENTIAL 


MIDDLE  OF 
INTERVAL 
0. 

1. 

2. 

3. 

4. 

5. 

1: 


NUMBER  OF 
OBSERVATIONS 

g  ********* 
10  ********** 

2  ♦* 

4  **** 


HISTOGRAM  NORMAL  (log  is  taken  then  histograa) 


MIDDLE  OF 
INTERVAL 


NUMBER  OF 
OBSERVATIONS 


1.0 

1 

* 

1.5 

3 

*** 

2.0 

6 

****** 

2.5 

6 

****** 

3.0 

5 

***** 

3.5 

2 

** 

4.0 

1 

* 

4.5 

2 

** 

-Figure  4  -  Histograas- 


B.  STATISTICAL  DISTRIBUTIONS 

i.  iks  SiESasniial  Pialsjjmiaa 

The  probability  density  function  of  the  exponential 
distribution  is  [Ref.  2], 

f(x,X)*Xe  ,  x  >  0  and  A>  0 

Its  cumulative  distribution  function  is, 

^dt-l-e 

which  can  be  easily  evaluated. 

The  exponential  distribution  is  used  frequently  as  a 
tine  to  failure  aodel  for  a  systea  when  a  constant  failure 
rate  is  assuaed  (Ref.  4].  But  for  repair  tiaes,  it  can  be 
shown  that  it  is  not  an  appropriate  aodel.  Because  repair 
tiae  includes  diagnostic,  correction,  and  verification  tasks 
a  repair  time  density  au3t  have  a  value  of  zero  at  tiae  t*0. 


It  should  then  increase  to  its  aaxiaua  value  rapidly  and 
then  gradually  decrease  towards  zero  as  tiae  increases. 


The  maximum  likelihood  estimator  of  the  parameter  of 
the  exponential  distribution  from  givan  data  is. 


X.-!- 


where 


X.  *U> 

Is  I 


2.  ifiafiaiaal  ^ion 

This  distribution  has  many  different  shapes  for  non¬ 
negative  variates.  It  is  skewed  to  the  right,  the  degree  of 
skewness  increasing  with  increasing  values  of  standard  devi¬ 
ation.  The  mean  and  standard  deviation  are  scale  and  shape 
parameters  respectively. 

A  random  variable  is  said  to  have  a  lognormal  dis  - 
tribution  if  the  logarithm  of  the  variable  is  normally  dis¬ 


tributed. 


XCJ27T 


with  mean  and  variance 
fl  *  B  (lnX)  an  d 


"  '  (Inx-u'f 

e  2  \  G  )  dx 


CT  *  Variance  (lnX) 


If  data  X(1),  X  (2)  ,  .  ,  X(n)  and  if 

T  ( 1 )  *lnX(1)  ,  T  (2)  =1  nX  (2)  ,  ....,  Y(n)*lnX(n)  than, 
maximum  likelihood  estimators  of  the  lognormal  distribution 


A 

11*1 


A  2  1 


(7  (-) -Y) 


The  p 

percentile  is 

calculated 

Hedian 

Z(.5)*0 

r,.5 

th 

90  per. 

Z(.9)»1.282 

x 

0.9 

th 

95  per. 

Z  (.9  5)  *  1.  6  45 

X 

0.95 

/1  +  1.  282  G* 

fj  +  1.645^ 


3 .  2a JIM  Distribution 

I  random  variable  x  has  a  gaama  distribution,  and  1 
is  referred  to  as  a  gamma  random  variable,  if  and  only  if 
its  probability  density  is  given  by,  [Ref.  9] 


f(*) 


b  n&) 
0 


a-1  -x/b 
-x  e 


for  x  >  0 


elsewhere 


where  a  2  0  and  b  >  0 

a  is  shape  parameter 
b  is  scale  parameter 
f^is  the  gamma  function 

The  mean  and  variance  of  the  gamma  distribution  are 
given  by, 
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aean=a.b  and 


variances,  b 


Special  cases  of  the  gamma  distribution  play  an  im¬ 
portant  roles  in  statistics.  For  instance,  for  a*1  and  b*Q 
we  get  the  exponential  distribution  with 


For  the  gaaaa  distribution  the  scale  and  shape  para- 
aeters  can  be  estiaated  fro  a  data  in  many  ways;  here  are  two 

a.  Method  Of  Moment  (H.O.N)  estimators 

2 

(sample  mean)  sample  variance 

a=— - .  b  - - - - . . . 

sample  variance  sample  mean 

2 

where  aean  =  a.b  and  variance=a. b 

b.  Mazimua  Likelihood  Method  (M.L.E.)  estimators 
This  aethod  is  more  complex  as  compared  to  the 

M.O.H.;  it  requires  nuaerical  calculation.  There  is  an  APL 
function  naaed  ESTGAM  that  gives  directly  M.L.E.  (App.C). 

It  is  used  because  it  provides  a  batter  estimator  of  “he  pa¬ 
rameters  a  and  b  in  large  samples,  provided  the  data  truely 
cone  from  a  gaaaa  model. 

C.  PBOBABILI TY  PLOTTING 

Plotting  data  points  on  probability  paper  is  quite  easy 
and  does  not  require  aany  calculations.  It  is  a  visual  exa¬ 
mination  rather  than  statistical  calculation  (Figure  5) . 

What  can  we  learn  from  a  probability  plot? 

1.  Visual  indication  of  assumed  distribution. 

2.  Examine  departures  froa  assumed  distribution. 


Ea.  Tatis  longer  than  normal 
b.  Tails  shorter  than  normal 
-Figure  5- 

i.  iesaoiaaJL  aisliifesitioa 

Consider  the  order  statistics  of  a  sample 

x(1)  <  x  (2)  <  . <  x  (N)  Then,  x 

is  an  estimate  of  the  (i/N)  th  percentile  of  X 
If  X  is  H(0,i)  then  we  expect 
-1 

x(i  )2zJf  (i/N)  for  1  =  1, 2, 3 . N 

Then  a  plot  of  X  (i)  vs  Jt  1  (i/H)  should  give  a 
straight  line.  The  only  problem  is 

-1  -1 
St  fi/fW  1 “infinity,  so  we  use  St  (i/(S+s)),  s  =  1. 
PROCEDURES 

a.  Use  normal  probability  paper 
Plot  (i/N+s)  ws  X(i> 

b.  Use  r.ormplot  in  APL  workspace  [Ref.  21] 
Use  MINITAB 

_ N SCORES  Of  Cl, put  in  C2 

PLOT  C2  VS  Cl 


NSCORES  are  (i/N+s> 


Consider  log  of  data  1(1)  <  Y  (2)  < 

. S  Y  (N)  where  Y(i)=LnX(i)  and  then  apply  above 

procedure  to  get  lognormal  plot. 

2.  ii£c aeaiial  eislsibution 

Suppose  X  is  exponential  A  with  sample  size  n,  so 
Pr(X  <  x) *  1-exp  (—A X)  where  X  >  0.  Then  (i/n+s)  is 

an  estimate  of  Pr  (X  <  X  (i) )  =1-exp  (-Ax  (i) ) 

So  if  the  data  really  cone  from  an  exponential  dis¬ 
tribution  (i/H+s)  vs.  1-e xp  (-Ax  (i) )  ,  should  yield  a 
straight  line. 


PROCEDURES 

a.  In  (i/H+s)  vs.  X  (i)  oi  rectangular  paper, 

b.  EXPOHPLOT  in  APL  workspace,  [Ref.  21] 

C.  MIHITAB, 


d.  Use  Gamma  probability  paper  for  ALPHA*1 
BET A*Q  where  Q  is  exponential  parameter  so  that, 

f  (x)*  (1/Q)  *exp(-X/3)  for  X  >  0 
3.  Gamma  Distribution 

There  is  an  APL  program  function  named  3A8  that  uses 


estimated  gamma  parameters  and  gives  test  statistic  D.  The 
program  also  gives  a  matrix  for  probability  plotting  pur¬ 
poses.  (see  Appendix  C) 


D.  STATISTICAL  TESTS  FOR  D ISTR IBOTIDRAL  ASSUMPTIONS 

He  need  to  know  whether  we  are  actually  dealing  with  a 
saaple  fro a  an  exponential  population  or  whether  the  data 
values  of  randoa  variables  have  a  logaoraal  or  3oae  other 
distribution.  A  statistical  test  of  significance  of  a  dis¬ 
tributional  assuaption  provides  an  objective  technique  for 
tentatively  assessing  whether  or  not  an  assuaed  aodel  pro¬ 
vides  an  adequate  description  of  the  observed  data. 

There  are  three  basic  steps  involved  in  statistical  test 
aethods. 

1.  A  test  statistic  is  calculated  from  the  data, 

2.  The  probability  of  obtaining  the  calculated  test 
statistic  is  deterained, 

3.  Assessaent  is  aade  of  the  adequacy  of  the  assuaed 
distribution. 

a.  If  the  probability  of  obtaining  the  calculated 
test  statistic  is  "LOW",  one  can  say  that  the  assuaed  dis  - 
tribution  does  not  provide  ar.  adequate  representation. 

b.  If  the  probability  of  obtaining  the  calculated 
test  statistic  "HIGH”  ,  then  the  data  provides  no  evidence 
that  the  assuaed  distribution  is  not  adequate,  judged  by  the 


The  definition  of  "LOH"  or  "HIGH"  depends  on  the 
user’s  preferences,  and  the  consequences  of  rejecting  the 
distribution.  For  example,  the  rejection  statistic  is  typi¬ 
cally  be  0.01,  0.05,  or  0.1.  For  the  purpose  of  this  study 
0.05  was  selected  as  the  reject  criterion. 

It  should  be  pointed  out  that  a  statistical  test, 
although  it  allows  one  to  reject  an  assumption  as  inadequate 
does  not  prove  that  the  assumed  distribution  is  correct. 

1 •  The  Chi-Squared  Soodness-of-?it  Test 

This  test  is  one  of  the  oldest  and  most  commonly 
used  for  evaluating  distributional  assumptions.  Basically, 
the  given  data  are  grouped  into  cells  and  compared  to  the 
expected  number  of  observation  in  the  cells  based  on  the 
assumed  distribution.  Then  if  the  test  statistic,  calcula¬ 
ted  from  this  comparison  exceeds  a  critical  chi-square  value 
the  assumed  distribution  is  rejected. 

The  only  problem  is  dividing  the  data  into  cells. 

If  the  number  of  observations  is  small  it  is  usually  sugges¬ 
ted  to  use  a  number  of  cells  as  large  as  possible,  subject 
to  the  restriction  that  it  should  not  exceed  N/5  (  M  is  sam- 


The  computations  involved  in  the  Chi-Squared  Good  - 
ness-Of-Fit  test  were  mads  by  a  computer  program,  in  which 
are  given  as  inputs  the  assumed  cumulative  distribution 
function,  the  number  of  observations,  the  number  of  equipro 
bable  cells,  and  the  number  of  parameters  estimated  from  th 
sample.  The  outputs  are  the  Chi-Square  statistic  and  its 
probability  of  exceeding  a  chi-square  value  for  a  given  num 
ber  of  degrees  of  freedom. 

The  procedure  used  is  described  as  follows:' Ref .  2] 
a.  The  cell  boundaries  are  determined  from  the 
assumed  cumulative  distribution  as  the  values  such  that  the 
estimated  probability  of  the  observation  value  falling 
within  a  given  class  is  1/K.  Por  each  class, 

Pr(X<x(i))  =  i/k,  where 

x*the  random  observation  to  be  assigned  to  the 

i-th  cell 

x  (i)  =the  i-th  cell  boundary  to  be  solved  from 
the  above  formula 

k  *  the  number  of  cells,  in  this  study  k  =  n/5 

The  lower  bound  of  the  first  cell  and  the  upper 
bound  of  the  last  cell  are  the  smallast  and  the  largest 
values  of  the  observations. 


b.  The  expected  number  of  observations  for  each 
cell  E(i)  is  equal  to  n/k  for  each  cell. 

c.  The  number  of  observed  values  in  each  cell  S(i) 
is  counted  based  on  the  results  of  the  above  equation. 

d.  The  test  statistic  is 


k  T  k  2  T 

*~rl  !i a<i*  l 


n 


2 

e.  The  computed  value  \  is  used  to  compute  the 
level  of  significance,  or  the  probability  of  a  Chi-Square 
value  with  degrees  of  freedom  equal  to  K-H-1  where  M  is  the 
number  of  estimated  parameters  from  the  sample  to  exceed 
calculated  chi-square 


2  2 

o(»Pr(  (V)  >  X  )  with  V  degrees  of  freedom 

1  -* 

If  oC  is  less  than  or  equal  to  .05  the  assumed  dis 
tribution  can  be  rejected  as  inadequate. 


2.  The  Kolmoqoro v-Smirnov  Test  ( K- S  Test) 

K-S  test  examines  the  cumulative  distribution  func¬ 
tion  rather  than  probability  density  function  .  In  simply 
compares  the  sample  cumulative  distribution  function  S(.) 
to  the  hypothesized  distribution  function  F(.).  As  a  mea- 


26 


sure  of  comparison  the  test  uses,  D  *  Max |  F  (. ) -S  (.)  |  ;  this 


is  the  main  difference  between  the  Chi-Squared  3.0.?.  test 
and  the  K-S.  test.  [Bef.  4] 

On  the  other  hand  here  we  are  dealing  with  Type  I 
error,  that  is  the  hypothesis  is  true  and  reject  it.  The 
probability  of  commit ing  a  Type  I  error  is  called  the  signi¬ 
ficance  level  for  the  test  alpha.  Here  alpha  is  taken  to  be 
0.05.  This  means  no  more  than  a  5  %  chance  of  rejecting  rhe 
true  hyphothesis. 

X  has  CDF  F  (x)  =  Pr  ( X  <  x)  , 

Take  a  sample  X(1)  <  X(2)  < . <  X(N), 

Be  want  to  estimate  F(x)  =  Pr(X  <  x)=  S  (x) 

n 


S  (x) 
n 


for  X  <  X(1) 


n 

1/N 

"  X 

> 

X  (1)  .  AMD.  X 

< 

xm 

VI 

If 

2/M 

"  x 

> 

X  (2)  .  AMD.  X 

< 

X  (35 

II 

• 

• 

• 

•  •  •  •  •  • 

•  •  •  • 

19 

ff 

K/M 

'•  X 

> 

X  (K)  .  AMD.  X 

< 

X(K+1) 

tl 

ft 

• 

• 

1 

• 

• 

X 

• 

• 

> 

i(M) 

•  •  •  • 

•  •  •  • 

lev  el , 


Find  D*  MAX|  S(x)  -F  (x)  | 

Using  Table  9  of  reference  4  for  given  significance 

If  D  exceeds  the  tabulated  value  that  comes  from  the 


table,  reject  the  hypothetical  distribution 


If  D  does  not  exceed  the  tabulated  value  that  comas 
from  table,  do  not  re  jeer  the  hypothetical  distribution. 

There  are  two  kind3  of  tables  for  the  K- 3  test.  One 
of  these  is  used  if  the  parameters  are  estimated  from  data, 
and  the  other  one  is  used  if  the  parameters  are  given.  Here 
parameters  are  estimated  from  data. 

3.  2fee  S  test  for  L  ogq  orma  1  Assumption 

The  3  test  is  shown  in  Ref.  3  to  be  an  effective 
procedure  for  evaluating  the  assumption  of  normality  against 
nonnormal  alternatives,  even  if  only  a  relatively  small  num¬ 
ber  of  observations  are  available.  Hahn  and  Shapiro  suggest 
that  the  W  test  may  also  be  used  to  evaluate  the  assumptions 
of  the  lognormal  distribution.  This  follows  because  of  the 
property  that  if  the  logarithms  of  the  observations  follow  a 
normal  distribution,  then  the  original  values  of  the  obser  - 
vations  are  lognormally  distributed. 

The  following  procedure  was  applied  using  tables 
from  Reference  2. 

a.  The  observations  are  ordered, 

X  (1)  S  X  ( 2)  S  ...<  X(i) . <  X(n-1>  <  X  (r.) 


b.  The  following  parameters  are  computed. 


1. 


2  p  ,2 

S  =  r  [lnx(i|] 


.X  lnx(i) 

i=1 


ii.  If  n  is  even,  k=n/2 

n  is  odd,  k=  (n- 1)  /2 

B*  .^C  a  (k-i+ 1)  M[  (lax  (n-iel) )  -  (lnx  (i) )  ]  ] 
where  the  values  of  a(k-i  +  1)  for  i=1 ,2, . .k 

are  given  in  table  IX  [Ref.  2]  for  n=3  to  r.*50 

iii.  The  test  statistic  W  is 

2  2 
»*  B  /  S 

c.  The  approximate  probability  of  obtaining  the 


calculated  value  of  W  can  be  obtained  from  the  formula 


v  ,  w-gl 

Z=A+7iln  - 1 

"  1-w( 

using  the  values  of  A  ,7^  and  £  given  in  Table  XI  of  Ref.  2 
for  the  appropriate  sample  size  and  then  using  standardized 
normal  distribution  to  determine  the  probability  of  obtain¬ 
ing  a  value  less  than  or  equal  to  Z,  which  is  the  signifi¬ 
cance  level  of  the  test. 


q  *  Pr(Z  S  Z) 

<3 

If  q  is  less  than  or  equal  to  .05,  the  selected  le¬ 
vel  of  significance  in  this  study,  the  lognormal  is  rejected 


as  an  inadequate  assumption 


III. DATA  ANALYSIS  BESDLTS 
A.  BESOLT  OF  TESTS  FOB  DIS TBIBUTIONAL  ASSUMPTIONS 

i.  sajaaasz  al  cssalts 

Tables  2,3  and  4  summarize  the  results  of  the  sta  - 
tistical  test  analysis.  In  all  cases  the  5  %  level  vas  cho¬ 
sen  as  the  reject  criterion.  The  46  sets  cf  data  used  in 
this  study  can  be  divided  into  three  parts  as  indicated  in 
Section  on  IB2. 

The  following  is  an  example  of  the  computer  program 
summary  table  which  includes  the  results  for  the  Chi-Square 
test  and  calculated  parameters  from  the  sample  data. 

SET  NO  18  -  AN/GXS-2  (?) 


SAMPLE  SIZE 

N  *  26 

NO.  OF  CELLS 

K  »  5  (a) 

SAMPLE  MEAN 

*  19.70 

STANDABD  DEV 

*  22.62  (b) 

EXPONENTIAL 

LOGNOBMAL 

DIFFEBENCE 

PABAH1 

0.05 

2.56 

<c) 

PABAM2 

0.78 

MTTB 

19.70 

19.  13 

2.96  *  (d) 

50 -TH  PEBCNT 

13.65 

12.93 

5.61  % 

90-TH  PEBCNT 

45.35 

40.22 

12.77  % 

95-TH  PEBCNT 

59.00 

55.46 

6.39  * 

CHI-SQB  ST  AT 

5.9  2 

0.54 

(e) 

DEG  OF  FBEED 

3 

2 

(f) 

SIGNIF  LEVEL 

0.  115  E+00 

0 • 764B+00 

<g> 

0 


zafel f  3 


LOGNORHAL  ABO  EXPONENTIAL  RESULTS 


I'rPT'TT 

P(X2)  | 

00166 


Q.48E-4 

0.59E-7 

0.018 

0.245 

0.238E-6 

0.0101 


ITTTm 
P(X  )  I  P(») 

"U7B59-1 - or 

0.358  0. 

0.05  88  0. 

0.0005  0. 

0.0128  0. 

0.0468  0. 

0.596E-7  — 

0.0143  0. 

0.459  0. 


K-S  Test. 
Ts.St.  Tab. 7a. 


T3S 

025 

0745 

0 

0328 

0002 


0.017 

0.31 


10 

50 

0.5  96E-7 

0.0445 

0.099 

11 

50 

0.00011 

0.979 

0.632 

12 

43 

0.596  E-6 

0.0421 

0.19 

13 

30 

0.28E-5 

0.01  86 

0.01 

14 

26 

0.  129E-4 

0.292 

0.238 

15 

50 

0.208 E-4 

0.00408 

0.001 

15a 

33 

0.2  16  E-4 

0.  484 

0.06 

16 

30 

0.8 11 E-4 

0.0321 

0.0 

17 

50 

0.0323 

0.494 

0.583 

18 

26 

0.  115 

0.764 

0.36 

19 

57 

0.596E-7 

0.360 

20 

57 

0.0 

0.  11 6  E-4 

21 

157 

0.0 

0.0 

22 

28 

0.381 

0.526 

0.624 

23 

36 

0.808 E-4 

0.  173E-3 

0.0074 

24 

58 

0.132E-3 

0.0255 

25 

20 

0.0334 

0.91 2E-2 

0.35 

26 

92 

0.6  56  E-6 

0.380  E-4 

27 

21 

0.  174E-2 

0.231E-2 

0.  175 

28 

29 

0.018  ~ 
0.69E-2 

8:81 

8:  III 

30 

35 

0.254 E-4 

0.592 

0.617 

31 

81 

0.186E-3 

0. 86 1 E-4 

32 

22 

0.468E-2 

0.06  3 

0.59 

33 

45 

0.4  27E-3 

0.315E-3 

0.013 

34 

171 

0.596E-7 

0.596E-7 

35 

31 

0.456 

0.724 

0.97 

36 

37 

0.0258 

0.571 

0.11 

37 

60 

0.542Z-5 

0.430E-2 

38 

20 

0.0334 

0.094 

0.896 

39 

95 

0.298  E-6 

0.22 

40 

22 

0.0716 

0.201 

0.424 

41 

35 

0.0476 

0.126 

0.48 

42 

51 

0.  193 

0.  103 

43 

38 

0.4  70E-2 

0.242 

0.04 

44 

87 

0.182E-2 

5.193 

45 

22 

0.468E-2 

5.340 

0.57 

46 

47 

0.7  57  E-5 

0.192E-3 

0.004 

TTUO’ 

.079 

All 

.148 

.27 

.172 

.235 

.217 

.212 

.086 

.143 

.3 

.151 
.245 
.211 
.199 
.125 
.084 
.23 
.441 
.073 
.145 
.132 
.143 
.16 
.102 
.301 
.143 
.106 
.263 
.121 
.219 
.147 
.069 
.117 
.083 
.123 
.099 
.091 
.156 
.  125 
.092 
.265 
.106 
.224 
.179 


TT35- 
.135 
.161 
.135 
.154 
.  193 
.115 
.190 

:??! 
.125 
.135 
.161 
.174 
.125 
.154 
.161 
.125 
.174 
.117 
.117 
.07 
.167 
.  147 
.116 
.  190 
.092 
.  193 
.154 
.130 
.149 
.098 
.189 
.132 
.068 
.159 
.145 
.114 
.110 
.091 
.189 
.149 
.124 
.  144 
.095 
.ie9 
.13 


notes: 


a.  The  number  of  e quiprobable  cells,  K,  was  chosen 
as  H/5  where  N  the  saaple  size. 

b.  The  saaple  Bean  and  standard  deviation  are  cal¬ 
culated  based  on  the  maximum  likelihood  estimates. 

c.  PARAM 1  for  the  exponential  distribution  is  the 
reciprocal  of  the  saaple  lean,  for  the  lognormal  distribu¬ 
tion  PAP  AMI  and  PAS  AM  2  are  /}  and  G ^  the  parameters  of  the 
lognoraal  distribution. 

d.  The  HTTP  and  the  50th,  90th  and  95th  percentiles 
are  calculated  froa  the  saaple  and  ace  based  on  the  rela¬ 
tionships  between  the  calculated  parameters  and  their  dis¬ 
tribution  functions.  The  percentage  difference  is  between 
the  exponential  and  lognormal  HTTP  and  percentiles. 

e.  The  Chi-Square  statistic. 

f.  The  number  of  degrees  of  freedom  (K-N-1)  is  K-2 
for  the  exponential  distribution  and  K-3  for  the  lognoraal. 
This  is  because  one  parameter  is  estimated  from  the  sample 
in  the  exponential  case  and  two  in  the  lognormal  case. 

g.  The  level  of  significance  is  the  probability  of 
a  Chi-Square  variate  with  the  specified  degrees  of  freedom 
exceeding  the  calculated  Chi-square  statistic. 


For  the  exponential  aodel,  only  the  Chi-Squared  test 
was  used.  In  6  out  of  46  cases  we  would  accept  the  exponen¬ 
tial  aodel.  If  we  coapare  the  exponential  vs.  lognorial  ao¬ 
del  using  the  Chi-Square  test,  test  prefers  the  lognoraal 
aodel  tc  the  exponential  aodel.  In  only  one  case  would  the 
exponential  aodel  be  accepted  and  the  lognoraal  rejected. 


lognoraal 


A 

S 


exponential 
1  B 


L5  . 

JLI 

i 

21  1 

6  40 


24 

22 


Therefore  the  lognoraal  aodel  is  a  better  descriptor 
with  respect  to  the  exponential  aodel. 

In  ao st  cases  the  lognoraal  assuaption  works  reason¬ 
ably  well.  In  those  cases  in  which  the  results  of  the  chi- 
squared  test,  W  test  and  K-S  test  do  not  agree,  a  plot  on 
lognoraal  probability  paper  was  used  to  determine  the  appro¬ 
priateness  of  this  assuaption  (*  in  Table  2).  Comparison 
between  test  methods  were  Bade  in  three  ways  for  the  lognor¬ 
aal  model  and  fabled  below. 


Chi-Sq 


A 

a 


18  28 


25 

21 


b.  Chi-Sq.  vs.  H  test  (34  cases) 


w 

A  R 


Chi-Sq. 


1  19 

-Li 

1  T 

9  1 

23 

1  1 

-S  test  (34 

K-S 

A  R 

L  2L 

JL1 

1  3 

8  1 

16  1  8 

For  the  lognormal  aethods  the  test  methods  sometimes 
give  different  results.  The  K-S  test  and  the  chi-square 
test  agree  in  74  %  of  the  cases.  The  chi-square  agrees  with 
the  H  test  in  82  S  of  the  cases.  R  and  K-S  test  agree  62  %. 
Therefore,  there  is  no  significant  difference  among  the 
three  test  aethods. 

A  previous  study  has  been  done  using  24  data  sets  of 
electronic  systeas  and  equipments,  U3ing  only  the  R  and  chi- 
squared  tests  for  the  lognormal  model  [Ref.  1].  The  table 
below  gives  a  comparison  between  the  exponential  and  lognor¬ 
mal  for  the  previous  study.  For  the  exponential  model  do 
not  reject  (A)  or  reject  (R)  decision  was  made  cased  on  the 
chi-squared  test.  For  the  lognormal  model  accept  by  either 
the  chi-square  or  R  test  was  used  as  a  criterion. 


lognor  sal 


A 

B 


exponential 
A  R 


21 

3 


8  16 

The  tables  below  show  same  kind  of  results  for  the 
current  study.  Only  difference  between  the  two  studies  is 
there  are  different  kind  of  systems  and  equipments  (e.g. 


mechanic  or  electronic)  in  the  current  study. 


ex  ponential 


A  |  5  i  23  |  28 

I0’”01"*1  B  lj“lj!li  18 

6  4  0~ 

Por  all  cases  (46  cases) 
exponential 


lognormal 


A 

LI 

I  7  1 

8 

R 

I  o 

-nrl 

5 

i 

1  2 

For 

electronic 

equipments  only 

(13  cases) 


These  results  indicate  that  the  lognormal  model  is 
better  than  the  exponential  model  for  repair  time  data. 

For  the  gamma  model,  the  HLE  estimator  is  more  po  - 
werful  than  the  MOB  estimator,  with  the  K-S  test  used  for 
testing  the  data.  The  decision  was  made  to  accept  the  re  - 


suits  of  the  SLE  when  the  two  estimators  gave  different 


results,  since  in  no  case  did  the  BOM  suggest  acceptance 
when  the  MLB  said  to  reject.  In  the  eight  cases  where  there 
was  a  difference,  therefore,  the  BLE  was  used  as  the  crite  - 
rion.  This  gives  the  following  results  for  the  gamma  model. 

SOS 


A 

8 

A 

1  19 

1  8  1 

27 

MLB 

8 

1  o 

1  19  | 

19 

19 

27 

Therefore 

using 

the  MLB, 

ws 

in  27  out  of 

46  cases .  In 

16 

lognormal  model  is  also  accepted  and  in  5  additional  cases 
all  three  distribution  assumptions  would  be  accepted.  Again 
only  one  case  (case  6)  would  we  accept  the  gamma  and  expo¬ 
nential  model  but  reject  the  lognormal. 

We  can  also  compare  lognormal  and  gamma  model  using 
the  K-S  test. 


gamma 
A  8 


lognormal 


H 


LI3 

..Li 

18 

1  13 

15  1 

28 

26  2  0 

The  above  table  shows  that  the  K-S  rest  prefers  the 


gamma  mcdel.  27  out  of  46  cases  accepted  the  gamma  model, 
but  only  18  out  of  46  cases  accepted  the  lognormal  model. 


a 


Tables  6-9,  (Appendix  B)  summarize  the  results  for 
all  cases  and  for  each  of  the  three  different  rypes  of  data. 
Table  10  (Appendix  B)  gives  total  results  in  terms  of  sample 
sizes  of  20-29,  30-39,  40-50  and  over  50.  For  exponential 
assumption,  increasing  sample  siza  increases  rejection  per  - 
centages.  For  gamma  assumption,  increasing  sample  size  also 
increases  rejection.  In  lognormal,  increasing  the  sample 
size  decreases  rejection  up  to  sample  sizes  of  forthy,  then 
rejection  increases  rapidly. 

Because  of  the  large  number  of  data  sets,  only  some 
interesting  cases  of  the  data  analysis  and  probability  plots 
for  those  sets  are  discussed  in  the  next  section. 

2.  Discussion  of  speci flc  cases 

The  following  discussions  illustrate  some  cases 
which  are  of  special  interest.  The  data  analysis  results 
and  plots  are  given  in  Appendix  A. 

a.  Set  No.  18  -  AN/GXS-2(V) 

This  case  is  an  illustration  where  all  three 
models  are  accepted  (not  rejected)  for  all  tests.  Data  were 
collected  from  26  maintenance  tasks.  In  this  case,  as  in 
most  other  of  the  formal  maintainability  demonstration  (ca¬ 
ses  8-20) ,  the  tests  were  performed  in  accordance  with  the 


D.S.  military  standard  for  maintainability  demonstration 
(HIL-STD-471)  [Bef.  23]  in  which  sample  selection  was  ge  - 
nerated  from  knowledge  or  estimates  of  component  failure 
rates  and  tests  were  conducted  using  service  technicians  of 
the  appropriate  skill  levels  and  given  prior  training  on  the 
equipments. 

This  case  fits  the  lognormal,  exponential  and 
also  gamma  for  both  estimators.  All  three  tests  give  high 
significance  for  lognormali ty. 

b.  Set  No.  6  -  Reactor  Coolant  Pump 

This  case  is  the  only  one  the  exponential  model 
accepted  and  the  lognormal  model  rejected.  However  the 
plots  are  not  good.  Further  investigation  is  necessary  in 
this  case  but  since  only  the  data  was  available,  it  was  not 
possible  for  us  to  do  so. 

c.  Sets  No.  28  and  41  -  Hain  Rotor  Head  and  Blades 
These  two  cases  indicate  a  very  good  example  of 

lognormality.  All  three  tests  signal  "do  not  reject"  for 
the  lcgnormal  assumption.  The  plots  are  also  very  good. 

The  gamma  family  does  not  well  represent  the  data  sets. 

d.  Sets  No.  1  and  2  -  French  Fessenheim  Pumps 
(  Down  /  Repair  )  times 
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These  two  cases  are  discussed  as  unsuccessful 
for  gamma  and  exponential  but  good  for  lognormal  since  at 
least  2  out  of  3  tests  give  acceptance  for  the  lognormal 
assumption.  For  repair  times  the  9  test  rejects  the  assump¬ 
tion  while  the  Chi-squared  and  K-S  tests  do  not,  but  down 
times  are  rejected  by  only  the  K-S  test.  There  are  some 
outliers  around  5X,  but  the  rest  of  the  data  show  a  good 
lognormal  fit  in  the  probability  plot. 

The  existence  of  outliers  could  be  a  reason  for 
rejection  by  one  of  the  3  tests.  For  downtime  data,  after 
trimming  the  two  high  outliers  the  lognormal  assumption  was 
also  accepted  by  the  K-S  test.  For  repair  time,  after  trim¬ 
ming  one  outlier  the  data  shows  high  significance  (0.15)  for 
the  9  test.  Thus,  in  some  cases  it  would  be  valuable  to 
examine  in  more  detail  the  reason  for  the  high  or  low  outli¬ 
ers.  In  this  case,  only  the  data  was  made  available  so  that 
this  was  not  possible. 

e.  Set  Ho.  4  -  Heater  Drain  Pump 
Set  No.  5  -  Main  Feedwater  Pump 
Set  No.  4  illustrates  the  case  where  all  three 
tests  for  the  lognormal  model  are  rejected.  Here  outliers 
are  important  factors.  The  effect  of  outliers  can  be  seen 


easily  from  the  boxplot,  histograa  and  probability  plot.  No 
5  is  almost  the  same  as  no  4,  the  only  difference  between 
the  two  being  that  the  K-S  test  barely  accepts  the  lognor  - 
aality  in  set  no  5.  For  these  cases,  trimming  the  outliers 
did  not  change  anything  in  so  far  as  significance  levels  or 
acceptance  or  rejection  are  concerned. 

f.  Set  No.  9  -  IN/ ARC- 1 62  (Intermediate) 

Set  No.  10  -  AN/ARC-163  (Intermediate) 

These  two  cases  are  discussed  together,  because 
they  give  opposite  results.  Derailed  reports  exist  in  both 
cases.  Test  technicians  were  used  for  running  the  tests  and 
removing  the  faults.  [Refs.  12  and  13] 

The  exponential  assumption  is  rejected  in  both 
cases.  In  set  No  9  the  lognormal  assumption  is  rejected  only 
by  the  K-S  test  while  the  other  tests  give  high  signifi¬ 
cance.  The  probability  plot  and  boxplot  also  signal  that 
"do  not  reject"  is  appropriate  for  this  case.  The  gamma 
model  fits  this  set  for  both  estimators. 

In  Set  No.  10  only  the  W  test  accepts,  bun  the 
Chi-Squared  test  barely  rejects  (0.045)  and  the  K-S  test 
rejects  the  lognormal  model.  However  the  boxplot  and  proba¬ 
bility  plot  are  fairly  good  for  this  set.  An  examination  of 


p 


the  report  indicates  that  sone  bias  exists  by  the  use  of 
"average"  tiaes  for  soae  repair  tiae  segments  for  auch  of 
the  data  instead  of  normal  data.  Therefore ,  a  further  eza  - 
aination  of  this  data  would  be  in  order, 
g.  Set  No.  15  -  UN  /DTK-  14(C) 

Set  No.  15a  -  Revised  AN/UYK-14(C) 

Set  No.  15  illustrates  what  happens  when  the  da¬ 
ta  points  are  not  selected  at  random  and  thus  introduce  bi  - 
as.  In  this  case,  33  tests  were  selected  according  to  the 
sample  selection  procedure  in  HIL-STD-471A  [Hef.  23].  How¬ 
ever  in  order  to  have  50  data  points  according  to  the  test 
plan  17  additional  data  points  using  "average"  repair  tiaes 
for  several  components  from  a  previous  test  were  used.  This 
therefore  inxroduces  a  bias  in  the  data,  resulting  in  the 
initial  rejection  of  the  lognormal  model  (Table  3).  After 
taking  out  the  biased  data  points  the  analysis  shows  that 
the  lognormal  model  is  a  good  descriptor.  The  histogram 
does  not  look  very  good,  but  the  probability  does  for  the 
lognormal  model. 

This  example  then  shows  that  if  data  does  not 
fit  an  assumed  distribution  and  if  the  reason  is  bias,  then 
after  finding  a  good  reason  for  throwing  out  some  bias  do  - 
ints,  the  data  freguently  will  agree  with  a  lognormal  model. 


43 


h.  Sets  No  2 1  and  34  -  Fuselage  structure 


These  tvo  cases  are  example  of  large  sample  si  - 
zes  (157  and  171).  For  both  data  sets  all  distributional 
assumptions  were  rejected  by  all  test  nethods.  But  for  the 
lognormal  case  probability  plots,  boxplots  said  histograms 
are  fairly  good.  These  cases  should  be  further  invesriga  - 
ted. 


B.  EHR0B3  IN  CALCULATED  PARAMETERS 

i.  nisi  is  am  v&es  jssaaiag.  SiBgsgatjaj 

Eisitikatiaa  uaisii  s£  Laiaasaai 

The  percent  error  in  the  mean-time-to-rspair  (MTTR) 
is  calculated  as  follows. 


M 

log 

M 


a  -  m  i 
1 og  ezp  I 

E  *  . . . . 1 

n 

log 

a  the  lognormal  mean  s  2 
a  the  exponential  mean  * 


x  100 

<S  A  i 

^♦O.SG*^ 

1 


where. 


exp 

The  results  of  the  percentage  error  in  the  mean 
which  are  summarized  in  Table  5.  All  cases  haws  an  error 
less  than  18  X.  Therefore  since  the  inherent  availabilit 
formula  normally  used  is. 


y 


HTBF  1 

I  BTBP  ♦  MTTR  MTTR 

1  ♦ - 

MTBF 
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s 

lafels  5 

»•[ 

PERCENTAGE 

ERROR  IN  STIR 

8HEN  ASSUMING  AN 

EXPONENTIAL 

INSTEAD  OF  LOGNORHAL  DISTRIBUTION 

Set 

Sample 

Sample 

Error  (%) 

No 

He  an 

Bean 

|3-2| 

(EXP) 

(LOG) 

1  -1” 

L’  * 

-2- 

-3- 

3 

W’?.' 

1 

2 

119.51 
24.  0 

1:!! 

3 

121.63 

140.84 

13.63 

Lj  5 

57.  21 

48.6  9 

17.49 

1  9 

^  10 

7.23 

7.22 

0.  16 

9.35 

9.34 

0.08 

1 1 

43.05 

4  3.59 

1.23 

12 

15.75 

15.74 

0.02 

14 

15.73 

15.54 

1.25 

£  15 

15.  85 

15.59 

1.65 

K;  17 

17.42 

17.76 

1.93 

ii  18 

19.7 

19.1  3 

2.96 

3  22 

1.  87 

2.07 

9.73 

^  23 

1.55 

1.48 

4.9 

25 

1.  76 

1.78 

1.36 

27 

0.57 

0.3  1 

1.36 

28 

0.41 

0.28 

1.29 

29 

2.45 

2.44 

0.69 

1  1° 

1.7 

1.72 

1.45 

1.64 

1.66 

0.85 

I  35 

4.03 

4.57 

11.84 

?  36 

2.55 

2.39 

6.96 

38 

2.26 

2.3  3 

2.96 

39 

4.22 

4.1 

2.73 

4  0 

2.5 

2.57 

2.57 

4  1 

3.5 

3.5  1 

0.27 

1:'  4  2 

3.63 

3.6  6 

0.63 

*  43 

2.4 

2.4  4 

1.78 

|  44 

1.72 

1.87 

7.88 

V;  45 

2.  58 

2.60 

0.58 
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and  since  for  a  reasonable  availability  HTTB  «  MTBF  (e.g. 
HTTR  is  aeasured  in  hoars,  MTBP  is  100's  of  hoars),  then  it 
is  shown  that  the  exponential  assumption  for  HTTB  introduces 
negligible  error  in  availability. 

C.  CONCLUSIONS 

Froa  the  data  analysis  in  this  study  as  well  as  the  pre¬ 
vious  study  [Bef.  1],  it  is  concluded  that  the  lognoraal 
model  is  a  good  descriptor  for  repair  times.  28  out  of  46 
sets  show  that,  within  a  0.05  level  of  significance  for  at 
least  one  of  the  3  test  methods  ,  this  model  can  not  be  re¬ 
jected.  Similarly,  the  data  analysis  shows  that  the  assump- 
tion  of  the  exponential  distribution  should  be  rejected  in 
40  sets.  The  gaaaa  assumption  is  not  rejected  in  27  cases. 

The  percent  difference  in  the  HTTB,  when  assuming  expo¬ 
nential  distribution  when  the  true  distribution  is  lognoraal 
has  been  found  to  be  negligible  for  calculated  system  avai¬ 
lability. 

when  the  results  of  the  statistical  tests  indicate  oppo¬ 
site  conclusions,  probability  plots,  boxplots  and  histograms 
are  helpful  in  determining  an  acceptable  model  and  to  explo¬ 
re  anomalies  in  the  data.  Probability  plots  are  also  useful 
in  estimating  percentiles,  as  well  as  density  parameters. 
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SET  NO  1  -  FRENCH  FESSENHEI H  PUMPS  (DOHN  TIME) 


SAMPLE  SIZE 

N  *  43 

NO.  OF  CELLS 

K  =  7 

SAMPLE  MEAN 

*  119. 51 

STANDARD  DEV 

=  225.79 

EXPONENTIAL 

LOGNORMAL 

DIFFERENCE 

PARAM1 

0.0  1 

3.89 

PARA  M2 

1.66 

MTTR 

119.5  1 

112.44 

6.29  % 

50-TH  PERCNT 

82.84 

49.10 

68.70  % 

90-TH  PERCNT 

275.1  9 

255.75 

7.60  % 

95 -TH  PERCNT 

358.02 

408.08 

12.27  % 

CHI-SQR  STAT 

13.1  6 

2.42 

DEG  OF  FREED 

5 

4 

SIGNIF  LEVEL 

0.  219E-01 

3 . 659E+00 

BOX PLOT  SET1 

- * 

BOX PLOT  LNSET1 


0  0 


♦  I* 


47 


HISTOGRAM  SET1 


MIDDLE  OF 

NUMBER 

OF 

INTERVAL 

OBSSRV  ATI  ON  S 

10.00 

6 

*****  * 

15.00 

6 

*****  * 

20.00 

0 

25.00 

1 

* 

30.00 

3 

*** 

35.00 

2 

** 

40.00 

1 

* 

45.00 

0 

50.00 

0 

55.00 

4 

**** 

60.00 

4 

**** 

65.00 

0 

70.00 

0 

75.00 

0 

80.00 

2 

** 

HISTOGRAM  LNSET1 

MIDDLE  OF 
INTERVAL 


NUMBER  OF 
OBSERVATIONS 


1.5 

1 

* 

2.0 

0 

2.5 

11 

*********** 

3.0 

2 

** 

3.5 

6 

****** 

4.0 

8 

*****  *** 

4.5 

3 

*** 

5.0 

6 

****** 

5.5 

2 

** 

6.0 

2 

** 

6.5 

1 

* 

7.0 

1 

* 

8 
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SET  NO  2  -  PRENCH  FESSENHEIM  POMPS  (REP1IR  TIME) 


SAMPLE  SIZE 

N  =  43 

NO.  OF  CELLS 

ii 

SAMPLE  MEAN 

*  24.00 

STANDARD  DEV 

=  31.95 

EXPONENTIAL 

LOGNORMAL 

DIFFERENCE 

PARAM1 

0.0  4 

2.73 

PARAH2 

MTTR 

24.00 

0.73 

22.  17 

8.25  % 

50-TH  PERCNT 

16.64 

15.40 

7.99  % 

90-TH  PERCNT 

55.26 

46.00 

20.14  % 

95-TH  PERCNT 

71.90 

62.70 

14.67  % 

CHI-SQR  STAT 

19.35 

4.37 

DEG  OF  FREED 

5 

4 

SIGNIF  LEVEL 

0.  166E-02 

0 . 358E+00 

PLOT  SET2 

-I  ♦  I - 

*  O 

O 

0 

BOX  PLOT  LNSET2 


HISTOGRAM  SBT2 


NUMBER  OF 
OBSERVATIONS 

12  ************ 

24  ************************ 

3  *** 

1  * 

1  * 

0 

0 

1  * 

1  * 


HISTOGRAM  LNSET2 

MIDDLE  OF 

NUMBER 

OP 

INTERVAL 

OBSERV  ATI  ON S 

1.5 

4 

**** 

2.0 

8 

*****  *** 

2.5 

11 

*********** 

3.0 

11 

*********** 

3.5 

5 

***** 

4.0 

1 

* 

MIDDLE  OF 
INTERVAL 
0. 
20. 
40. 
60. 


80. 

100. 

120. 

140. 

160. 


SET  NO  4 

-  HEATBH  DRAIN  POMP 

SAMPLE  SIZE 

N  =  43 

NO.  OF  CELLS 

SAMPLE  MEAN 

=  67.60 

STANDARD  DEV 

EX  PONENTIAL 

LOSNORHAL 

PABAM1 

0.01 

2.90 

PAHA  M2 

1.81 

MTTB 

67.60 

44.72 

50 -TH  PERCNT 

46.86 

18.09 

90-TH  PERCNT 

155.67 

101.51 

95-TH  PERCNT 

202.53 

165.42 

CHI-SQR  ST AT 

56.47 

20.00 

DEG  OF  FBEED 

5 

4 

SIGNIF  LEVEL 

0.  596E-07 

0.5  00E-03 

BOXPLOT  SET4 

♦  10  0  0 

0  0  0 

K  *  7 

=  154. 16 
DIFFERENCE 


51.18  % 
>  100  % 
53.36  *5 
22.43  % 


BO t PLOT  INSET 4 


HISTOGEAH  SET4 


MIDDLE  OF  NUMBER  OF 
INTERVAL  OBSERVATIONS 

5.00  9  ********* 

10.00  12  ************ 

15.00  7  ******* 

20.00  6  ****** 

25.00  2  ** 

HISTOGRAM  LNSET2 

MIDDLE  OF  NOMBER  OP 
INTERVAL  OBSERVATIONS 


1.5 

5 

***** 

2.0 

9 

*****  **** 

2.5 

9 

********* 

3.0 

12 

************ 

3.5 

1 

* 

4.0 

1 

* 

4.5 

0 

5.0 

2 

** 

5.5 

0 

6.0 

2 

** 

6.5 

2 

** 

54 


SET  HO  5  -  MAIM  FEEDHAIEF  PUMP 


SAMPLE  SIZE 

N  =  33 

NO.  OF  CELLS 

SAMPLE  MEAN 

=  57.21 

STANDARD  DEV 

EXPONENTIAL 

LOGNORMAL 

PA5AM1 

0.0  2 

3.20 

PARAM2 

1.38 

MTTS 

57.21 

48.69 

50  -TH  PEBCNT 

39.6  6 

24.47 

90-TH  PERCNT 

131.74 

110.  10 

95-TH  PEBCNT 

171.39 

168.54 

CHI-SQB  ST AT 

11.91 

10.82 

DEG  OF  FREED 

4 

3 

SIGNIF  LEVEL 

0.  180  E- 01 

0. 128E-01 

BOX PLOT  SET5 

1*  I—  0  0  0 


K  =  6 

=  105.37 
DIFFERENCE 


17.49  ^ 
62.03  % 
19.65  *5 
1.69  * 


BOX PLOT  LNSBT5 


HISTOGRAM  SET 5 


MIDDLE  OF 
INTERVAL 
5.00 
10.00 
15.00 
20.00 
25.00 
30.00 
35.00 
40.00 
45.00 
50.00 
55.00 
60.00 


NUMBER  OF 
OBSEBV  ATIONS 
2  ** 

10  ********** 

3  *** 

3  *** 

0 

5  ***** 

1  * 

1  * 

2  ** 

1  * 


HISTOGBAM  LNSET5 


MIDDLE  OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

1.0 

1 

* 

1.5 

1 

* 

2.0 

2 

** 

2.5 

10 

********** 

3.0 

4 

**** 

3.5 

7 

*****  ** 

4.0 

4 

**** 

4.5 

0 

5.0 

0 

5.5 

2 

** 

6.0 

2 

** 

57 


SET  HO  6  -  REACTOR  COOLANT  POMP 


SAMPLE  SIZE 

N  *  21 

NO.  OF  CELLS 

SAMPLE  MEAN 

*  58.57 

STANDARD  DEV 

EXPONENTIAL 

LOGNORMAL 

PARA  Ml 

0.0  2 

3.64 

PARAM2 

1.49 

MTTR 

58.57 

80.08 

50-TH  PERCNT 

40.6  0 

38.05 

90 -TH  PERCNT 

134.87 

181.78 

95-TH  PERCNT 

175.46 

283.05 

CHI-SQR  ST AT 

2.81 

3.95 

DEG  OF  FREED 

2 

1 

SIGNIF  LEVEL 

0.  245  E+00 

0.468E-01 

BOXPLOT  SETS 


K  =  4 

*  39 . 22 

DIFFERENCE 


26.86  % 
6.69  % 
25.81  % 
38.01  % 


BOXPLOT  LNSET6 


HISTOGHAH  SET6 


SET  HO  9_- _AN/ARC-162  (7)  (INTERMEDIATE  LEVEL) 


SIMPLE  SIZE 

N  »  49 

NO.  OF  CELLS 

K  *  10 

SAMPLE  MEAN 

=  7.  23 

STANDARD  DEV 

*  3.44 

EXPONENTIAL 

LOSNORHAL 

DIFFERENCE 

PARAM1 

*  0.14 

1.38 

PAHA  M2 

0.19 

MTTR 

7.23 

7.22 

0.16  % 

50 -TH  PERCNT 

5.0  1 

6.57 

23.64  X 

90-TH  PEHCNT 

16.65 

11.49 

44.98  * 

95 -TH  PEHCNT 

21.67 

13.46 

60.99  X 

CHI-SQH  STAT 

47.  94 

6.71 

DEG  OF  FREED 

3 

7 

SIGNIF  LEVEL 

0.  119E-06 

3.4592+00 

BOX PL 07  SET9 

- 1  ♦  X -  ** 


BOX  PLOT  LNSET9 


HIS TOGBAH  SET9 


HIDDLE  OP 

NUMBER 

OP 

INTERVAL 

2. 

OBSERV ATI ON S 

0 

4. 

12 

************ 

6. 

18 

****************** 

8. 

6 

****** 

10. 

6 

*****  * 

12. 

3 

*** 

14. 

1 

* 

16. 

2 

** 

18. 

1 

* 

HISTOGRAM  LNSET9 


HIDDLE  OP 


NUMBER  OF 
OBSERV  ATIONS 


1.0 

2 

** 

1.2 

3 

*** 

1.4 

5 

***** 

1.6 

5 

***** 

1.8 

14 

*****  ********* 

2.0 

7 

******* 

2.2 

4 

**** 

2.4 

4 

**** 

2.6 

2 

** 

2.8 

3 

*** 
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SET  NO  10  -  AN/ARC-1 63  (7)  (INTERMEDIATE  LEVEL) 


SAMPLE  SIZE 

N  *  50 

NO.  OF  CELLS 

K  =  10 

SAMPLE  MEAN 

=  9.  35 

STANDARD  DEV 

=  4.27 

EX  PONENTIAL 

LOGNORMAL 

DIFFERENCE 

PARAM1 

0.1  1 

2.14 

PARAM2 

0.  18 

MTTR 

9.3  5 

9.34 

0.08  % 

50 -TH  PE  RC  NT 

6.4  8 

8.53 

24.03  % 

90-TH  PERCNT 

21.52 

14.73 

46.09  % 

95-TH  PERCNT 

28.0  0' 

17.20 

62.82  % 

CHI-SQR  STAT 

53.20 

14.40 

DEG  OP  FREED 

3 

7 

SIGNIF  LEVEL 

0.596E-07 

0.445E-O1 

BOX  PLOT  SET  10 


........  £  4  I- .... ........  4  44  4  44 

BOX  PLOT  LNSET10 


68 


HISTOGRAM  SET10 


I 

i 


MIDDLE  OF 
INTERVAL 
4. 

6. 

8. 

10. 

12. 

14. 

16. 

18. 

20. 


NUMBER  OF 
OBSERVATIONS 

5  ***** 

10  ********** 

15  *************** 

9  ********* 

2  ** 

1  * 

3  *** 

3  *** 

2  ** 


i 


HISTOGRAM  LNSET10 


MIDDLE  OF  NUMBER  OF 
INTERVAL  OBSERVATIONS 
1.2  1  * 


1.2 

1 

* 

1.4 

3 

*** 

1.6 

3 

*** 

1.8 

8 

*****  *** 

2.0 

7 

******* 

2.2 

15 

*****  ********** 

2.4 

4 

**** 

2.6 

1 

* 

2.8 

5 

***** 

3.0 

3 

*** 

SET  HO  15  -  AN/OYK-1 4  (C) 


SAMPLE  SIZE 

N  =  50 

NO.  OF  CELLS 

SAMPLE  MEAN 

»  15.85 

STANDARD  DEV 

EX  PON  ENT  I AL 

LOGNORMAL 

PARA  HI 

0.0  6 

2.58 

PARAH2 

0.33 

HTTR 

15.85 

15.59 

50- TH  PERCNT 

10.99 

13.21 

90 -TH  PERCHT 

36.50 

27.63 

95- TH  PERCNT 

47.4  8 

34.05 

CHI-SQH  ST AT 

35.60 

20.80 

DEG  OF  FREED 

8 

7 

SIGNIF  LEVEL 

0.  208  E-04 

0.408E-02 

BOXPLOT  SET15 

- 1  ♦ 

~ - 

—  * 

K  *  10 
=  11.06 
DIFFERENCE 


1.65  % 
16.85  X 
32.09  X 
39.45  X 


BOXPLOT  LHSET15 


SET  NO  15 A  -  REVISED  AN/JYK-14(C) 


SAMPLE  SIZE 

N  =  33 

NO.  OF  CELLS 

SAMPLE  MEAN 

*  14.37 

STANDARD  DEV 

EX  PONSNIIAL 

LOGNORMAL 

PABAH1 

0.0  7 

2.50 

PARA M2 

0.31 

MTTB 

14.37 

14.19 

50- TH  PEBCNT 

9.96 

12.15 

90 -TH  PEBCNT 

33.  10 

24.83 

95-TH  PEBCNT 

43.06 

30.40 

CHI-SQB  ST AT 

26.82 

2.45 

DEG  OF  FREED 

4 

3 

SIGNIF  LEVEL 

0.  216  E- 04 

0. 484E+00 

BOX  PLOT  SET  15  A 

i 

1  M 
1 

1 

1 

1  ♦ 

1 

1 

1  H 
1 

- -  * 

*  * 

K  =  6 

9.66 

DIFPEBENC3 


1.29  % 
17.98  * 
33.30  % 
41.65  % 


BOX PLOT  LN SET  15 A 


SET  NO  18  -  AN/GXS-2  (V> 


SAMPLE  SIZE 

N  =  26 

NO.  OF  CELLS 

SAMPLE  MEAN 

*  19.70 

STANDARD  DEV 

EXPONENTIAL 

LOGNORMAL 

PARAN1 

0.0  5 

2.56 

PARA  M2 

0.78 

MTTR 

19.70 

19.13 

50-TH  PERCNT 

13.65 

12.93 

90-TH  PERCNT 

45.35 

40.22 

95- TH  PERCNT 

59.00 

55.46 

CHI-SQS  ST AT 

5.92 

0.54 

DEG  OF  FREED 

3 

2 

SIGNIF  LEVEL 

0.  115B+00 

0.764E+00 

BOXPLOT  SET18 

— !*♦  I - 

O 

K  =  5 

*  22.62 
DIFPESENCE 


2.96  % 
5.61  % 
12.77  * 
6.39  r, 


BOX PLOT  INSET 18 


HISTOGRAM  SET18 

MIDDLE  OP 

NUMBER 

OF 

interval 

OBSERVATIONS 

0. 

2 

** 

10. 

14 

*****  ********* 

20. 

5 

***** 

30. 

2 

** 

40. 

0 

50. 

0 

60. 

1 

* 

70. 

0 

8  0. 

0 

90. 

2 

** 

HISTOGRAM  INSET  18 

MIDDLE  OP  NUMBER  OF 
INTERVAL  OBSERVATIONS 
1.0  1  * 

1.5  3  *** 

2.0  6  ****** 

2.5  6  ****** 

3.0  5  ***** 

3.5  2  ** 


Bet  V*o 
Giimm  P(«4 

0'  =  /.JJ 


ll 

<3. 


82 


SET  NO  21  -  FUSELAGE  STRUCTURE 


SAMPLE  SIZE 

N  =  157 

NO.  OF  CELLS 

SAMPLE  MEAN 

=  1. 54 

STANDARD  DEV 

EX  PON  ENTIAL 

LOGNORMAL 

PARA  Ml 

0.6  5 

-0.06 

PARAM2 

0.97 

MTTR 

1.54 

1.53 

50- TH  PERCNT 

1.07 

0.94 

90 -TH  PERCNT 

3.56 

3.33 

95 -TH  PERCNT 

4.6  3 

4.75 

CHI-SQR  ST  AT 

169.06 

155.0 

DEG  OF  PREED 

22 

21 

SIGNIF  LEVEL 

0.0 

0.0 

BOXPLOT  SET21 

— 1~~  1 - 

-  **  ** 

OOO  0  0 

K  =  24 

1.83 

DIFFERENCE 


1.01  % 
13.52  % 
6.93  % 
2.64  % 


0  0< 


BOX  PLOT  LNSET21 


HISTOGRAM  SET  21 

EACH  *  REPRESENTS 


2  OBSERVATIONS 


MIDDLE  OP 
INTERVAL 
0. 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 

10. 


NUMBER  OF 
OBSERVATIONS 


33 

73 

27 

7 

4 

2 

5 
3 
1 
0 
2 


*****  ****  ******  ** 

************************************* 

************** 

**** 

** 

* 

*** 

** 

* 


HISTOGRAM  LNSET21 


MIDDLE  OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

-2.5 

3 

*** 

-2.0 

2 

** 

-1.5 

12 

*****  ******* 

-1.0 

16 

**************** 

-0.5 

34 

*********  ************************* 

0.0 

36 

************************************ 

0.5 

27 

********* ****************** 

1.0 

10 

********** 

1.5 

7 

*****  ** 

2.0 

8 

*****  *** 

2.5 

2 

** 

SET  SO  28  -  MAIN  ROTOR  HEAD  AND  BLADES  (ELAPSE  TIME) 


SAMPLE  SIZE 

N  =  33 

NO.  OF  CELLS 

K  *  6 

SAMPLE  MEAN 

»  2.  07 

STANDARD  DEV 

=  1.87 

EX  PON  ENT  I AL 

LOGNORMAL 

DIFFERENCE 

PARAM1 

0.4  8 

0.46 

PABAM2 

0.52 

MTTR 

2.0  7 

2.05 

0.89  % 

50- TH  PERCNT 

1.44 

1.58 

9.09  f 

90 -TH  PERCNT 

4.77 

4.00 

19.31  % 

95 -TH  PERCNT 

6.21 

5.20 

19.33  * 

CHI-SQR  ST AT 

1 1.9  1 

2.09 

DEG  OF  FREED 

4 

3 

SIGNIF  LEVEL 

0.  180  2-01 

0. 554E+00 

BOXPLOT  SET28 

—I  ♦  1- 


BOXPLOT  LNSET28 


HISTOGRAM  SET 28 


HIDDLE  OF 
INTERVAL 
0. 

1. 

2. 

3. 

ft. 

5. 

6. 

7. 

8. 

9. 

10. 

11. 


NUMBER  OF 
OBSERVATIONS 
0 


16 

6 

7 

2 

1 

0 

0 

0 

0 

0 

1 


*****  *********** 
****** 

*****  ** 

** 

* 


HISTOGRAM  LNSET28 


HIDDLE  OF 

NUMBER 

OF 

INTERVAL 

OBSERVATIONS 

-0.8 

3 

*** 

-O.ft 

2 

** 

0.0 

8 

*****  *** 

o.ft 

6 

*****  * 

0.8 

6 

****** 

1.2 

6 

****** 

1.6 

1 

* 

SAMPLE  SIZE 

N  =  171 

NO.  OF  CELLS 

K  =  34 

SAMPLE  MEAN 

=  2.  83 

ST&NDABD  DEV 

*  4. 

60 

EX  PON  ENT  IAL 

LOGNOBMAL 

DIFFERENCE 

PABAM1 

0.3  5 

0.36 

• 

PAHA M2 

1.26 

MTTB 

2.8  3 

2.68 

5.26 

% 

50 -TH  PEBCNT 

1.96 

1.43 

36.77 

% 

90 -TH  PEBCNT 

6.5  1 

6.03 

7.95 

% 

95-TH  PEBCNT 

8.46 

9.05 

6.51 

% 

CHI-SQB  ST AT 

126.65 

92.85 

DEG  OF  FBEED 

32 

31 

SIGNIF  LEVEL 

0.  596  E-07 

0 . 5  96E-07 

BOXPLOT  SET34 


1+  i——*  **QO  0  000  0  0 


BOXPLOT  LNSET34 


HISTOGRAM  SET 34 


5  OBSERVATIONS  ABE  BELOtf  THE  FIRST  CLASS 


MIDDLE  OF  NUMBER  OF 
INTERVAL  OBSERVATIONS 


0.500 

38 

*****  ********************************* 

1.000 

33 

*****  **************************** 

1.500 

25 

*****  ******************** 

2.000 

15 

*****  ********** 

2.  500 

9 

********* 

3.000 

13 

************* 

3.500 

2 

** 

4.000 

3 

*** 

4.500 

3 

*** 

5.000 

3 

*** 

5.500 

0 

6.000 

5 

***** 

6.500 

1 

* 

7.000 

1 

* 

HISTOGRAM  LNSET34 


MIDDLE  OF 
INTERVAL 
-  2.  500 
-2.250 
-2.000 
-1.750 
-1.500 
-1.250 
-1.000 
-0.750 
-0.500 
-0.250 
0.0 
0.2 
0.5 
0.7 
1. 


NUMBER  OF 
OBSERVATIONS 
2  ** 

0 

0 

3  *** 

5  ***** 

0 

*********** 

15  *************** 

£  4i4i4tA  A  4i 

13  ************* 

19  ******************* 

15  *************** 

19  ******************* 

12  ************ 

IS  ****************** 

3  *** 

3  ***** *** 

5  ****** 

2  ** 

7  ***** ** 


SET  NO  41  - 

MAIN  ROTOR 

HEAD  AND  BLADES 

SAMPLE  SIZE 

N  =  35 

NO.  OP  CELLS 

SAMPLE  MEAN 

=  3.  50 

STANDARD  DEV  - 

EXPONENTIAL 

LOGNORMAL 

PABAM1 

0.2  9 

0.90 

PARAM2 

0.71 

MTTR 

3.5  0 

3.51 

50-TH  PEHCNT 

2.43 

2.47 

90 -TH  PEBCNT 

8.06 

7.25 

95-TH  PEBCNT 

10.49 

9.83 

CHI-SQR  ST AT 

11.20 

7.20 

DEG  OF  FREED 

5 

4 

SIGNIP  LEVEL 

0.476E-01 

0. 126E+00 

BOX PLOT  SET41 


I  ♦  I  — 


(MAH  HOURS) 

K  =  7 

3.42 

DIFPERENC 


0.27  % 
1.57  % 
11.29  % 
6.72  % 


BOX  PLOT  LNSET41 


i  w 


HISTOGRAM  SET41 


HI DDLS  0? 
INTERVAL 
0. 

2. 

4. 

6. 

8. 

10. 

12. 

14. 

16. 

18. 


NUMBER  OF 
OBSERVATIONS 
5  ***** 

13  ************* 

10  ********** 

4  **** 

1  * 

1  * 

8 

0 

1  * 


HISTOGRAM  LNSET41 


BIDDLE  OF  NUMBER  OF 
INTERVAL  OBSERVATIONS 


0.5 

3 

*** 

0.0 

7 

*****  ** 

0.5 

7 

******* 

1.0 

3 

*** 

1.5 

9 

********* 

5 

***** 

sstiaators  R 


-*J«  |  <N  I  00 


I able  2 

POMP  DATA  HESOLTS  (7  sets) 

Exponential  Distribution  :  Only  Chi  Square  test  used. 

Accept  Reject 


! 

I 

I 

j 

1 

1 

6 

1 

1 

Lognormal  Distribution  :  Chi  Square,  Kolmogorov  Smirnov 

and  (f  test  used. 


K-S 

test 

A 

w  tes 

4» 

•«. 

R 

W  test 

A 

R 

A  R 

K-s  test 

A  R 

1 

1 

1 

A 

Chi  Sq. 

1  o 

1 

1 

2  0 

Ch. Sq. 
test 

- 

test 

R 

1 

1 

0  2 

1 

Rl  0  |  1  1 

( 

6 

cases  ) 

(  1  case  ) 

Gamma  Dis 

tribution  : 

Kolmogorov 

Smirnov  test  used. 

M.O  . 

M.  a 

stimators 

A 

R 

M.L.  E. 

A  j 

0 

1  1  ! 

1 

i  i 

1 

estimators  R  I  0  j  6 


lafels  § 

BEPOBT  DATA  RESULTS  (13  sets) 

Exponential  Distribution  :  Only  Chi  Square  test  used. 

Accept  Reject 


1 

1 

12 

1 

_ l 

_ [ 

Lognormal  Distribution  :  Chi  Square,  Kolmogorov  Smirnov 

and  H  test  used. 


Chi  Sq. 
test 


K-S  test 


A 

W  test 
A  R 


R 

II  test 
A  R 

2  0 


0  0 


2  3 


(  11  cases  ) 


No  H  test 


K-S  test 
A  R 


Ch.Sq.  A  0 
test  — — 

R  0 


(  2  cases  ) 


Gamma  Distribution  ;  Kolmogorov  Smirnov  test  used. 


B.O.M.  estimators 
A  R 


H.L.E.  A 

estimators  R 


96 


l&lS  2 


FIELD  DATA  BE  SULTS  (26  sets) 


estimators  3 


labls  12 


TOTAL  X  RESULTS  (w.r. t.  sample  sizes) 


For  exponential  distribution  : 


Sample  Size 

A 

a 

20-29 

30.8 

59.2 

30-39 

8.3 

91.7 

40-50 

0. 

130. 

GT.50 

8.3 

91.7 

Fcr  lognormal  distribution  : 

sf?iie 

Chi- Sq.test 

a  a 

V 

A 

K-S  tes 
A  R 


20-29 

53.3 

46.2 

71.4 

28.6 

46.2 

53.8 

30-39 

66.7 

33.3 

58.3 

41.7 

50. 

50. 

40-50 

58.3 

41.7 

58.3 

41.7 

33.3 

66.7 

GT.50 

25. 

75. 

— 

— 

16.7 

83.3 

For  gamma  distribution  : 

Sample  A  in 


a 


in 


A  in  HLE 


*Emm  s 


COMPOTES  PBOGBAHS  FOB  DATA  ANALYSIS 

1.  LGHBML  (FORTRAN)  [Ref.  1] 

The  prograa  aakes  a  Chi -Squared  Goodness-Of-Fit  test  £or 
the  exponential  and  lognoraal  assuaptions.  The  aain  prograa 
computes  expected  values  of  ordered  observations,  saaple 
aean  and  standard  deviation  based  on  the  aaxiaua  likelihood 
estiaates  of  the  paraaeters  of  the  lognoraal  distribution, 
calculates  the  50th,  90th,  95th  percentiles  of  the  exponen¬ 
tial  and  lognoraal  distribution  assuaptions,  percentage 
differences  between  the  calculated  paraaeters  of  the 
exponential  and  lognoraal  distributions. 

2.  BOX  PLOT  (J5IN1TAB)  [Ref.  5] 

Syntax:  Boxplot  Vector 

Description:  Generates  a  boxplot  display  for  a  vector 
of  data.  A  rectangular  box  with  ends  corresponding  to  lower 
and  upper  quartiles  is  presented  with  the  aedian  marked  with 
an  plus  sign,  and  lines  are  drawn  on  each  side  of  the  box 
with  crosses  narking  the  lowest  and  highest  data  values  wit¬ 
hin  half  quartile  distance  in  each  side.  Outliers  are  nark¬ 
ed  with  stars  and  circles. 

3.  NORHPLOT  (APL)  [Ref.  21  ] 

Syntax:  Nornplot  vector 

Paraaeters:  Bid-  controls  the  horizontal  size  of  display 
Dep-  controls  the  vertical  size  of  display 
Subprograas:  BSCOBES,  PLOT 

Description:  Generates  a  normal  probability  plot  for  a 
vector  of  data.  If  data  fits  noraal  then  plot  is  straight 
line. 

4.  EXPONPLOT  (APL)  [Ref.  21] 

Syntax:  Exponplot  vector 

Paraaeters:  Bid-  controls  the  horizontal  size  of  display 
Dep-  controls  the  vertical  size  of  display 
Subprograms  PLOT 


Description:  Generates  a  exponential  probability  plot 
for  a  vector  of  data.  If  data  fits  exponential  then  plot 
is  straight  line. 

5.  G  AH  (A  PL) 

Syntax:  GAH  vector 

Sub  pro  grass:  E  ST  GAH,  GAHHA 

Description:  Per fora  K  olaogorov- Smirnov  Goodness-of-Fit 

test  and  gives  test  statistic  and  "reject"  /  "do  not  reject" 
idea  for  gaana  assuaption  for  both  HDR  and  HLE  estiaators. 

6.  KOL  (APL) 

Syntax:  Kol  vector 

Subprograa:  HOH,  IHVV 

Description:  Perforas  a  Kolmogorov  Smirnov  Goodnees-Of- 
Fit  test  for  a  data  vector  for  lognormal  assumption  and 
print  test  statistic.  If  test  statistic  is  bigger  than  the 
tabulated  value  "reject"  the  assumed  distribution,  otherwise 
"do  not  reject". 

7.  9  (APL) 

Syntax:  A  s  Vector 

Description:  Perforas  a  W  test  for  the  lognormal  as  - 
sumption.  A  must  be  a  data  vector  (from  Table  IX,  Hef.  2) 
and  X  vector  that  is  testing  against  lognornality.  A  and  X 
aust  have  the  same  dimension,  output  will  be  V  test  statis¬ 
tic.  Then  using  the  value  of  A, ft  and  6  (Table  XI,  Ref.  2) 
for  the  appropriate  sample  size  and  than  using  standardized 
normal  distribution  deteraine  the  probability  of  obtaining 
value  less  than  or  egual  Z,  which  is  the  significance  level 
of  the  test. 


9  O AM  C03? 

9  QAM  r 

C13  A  PROGRAM  PERFORMS  KOLMOGOROV  SMZRNOF  TEST  FOR 
C21  A  GAMMA  ASSUMPTION  FOR  MOM  AND  MLE  ESTIMATORS 

£33  i»,E*»,3 

[ 4  J  ESTGAM  T 

£53  Xlt-T-rECl  #23 

£63  X2+-Y+EC2J23 

£73  HffY 

£83  2«-o.886*h*o.5 

£93  AH*e[im 

£103  A2*,KC2J  13 

£113  »  (N-l)  )-rN 

£123  x«-u«-0 

£133  J<-JJ«-0/0 

£143  i~2:i«-i  +  1 

£153  -»li*ui>H) 

£163  A1  oamma  XIC13 
£173 

£183  -»L2 

£193  ‘-1JR«-|  <  J-K) 

£203  RR**£t*3 

£213  '  test  statistic  for  mom  estimators  is  *»fRR£i3 

£223  -»«-3M<RRm<2) 

£233  1  feject  for  alphA3<),05  for  mom  estimators  • 

£243  -H-5 

£253  «-3J  '  ®o  HOT  REJECT  FOR  ALPHAx0.05  for  mom  estimators 
£263  L5}ix«.ii  +  i 
£273  *»L6K\<xx>N) 

£283  A2  oamma  X2£ix3 
£293  JJf JJ ? f 
£303  -»*-5 

£313  <-6JS«-l  (JJ-K) 

£323  ss<.3C»s3 

£333  '  test  statistic  for  mle  estimator  is  1  rtssci3 

£343  *L7x\ (SSC13<2) 

£353  '  REJECT  FOR  ALPHAx0»05  for  mle  estimators  • 

£363  +99 

£373  L7j*  hot  reject  for  alphax0*05  for  mle  estimators 

£383  +99 

9 


GAM  SET 1 

TEST  STATISTIC  FOR  MOM  ESTIMATORS  IS  0.363 
REJECT  FOR  ALPHAs<).05  FOR  MOM  ESTIMATORS 
TEST  STATISTIC  FOR  MLE  ESTIMATOR  IS  0.154 
REJECT  FOR  ALPHAxO.05  FOP  MLE  ESTIMATORS 


Reproduced  from 
best  evoifobjo  copy. 


101 


cm 

m 

IM 

161 

Z71 

LB1 

ilh 

an 

C123 
C 133 
C 14] 
C 153 
C 16] 
C 17] 
C 18] 
C19] 
C20] 
C21  ] 
C22] 

m 

C25] 
C26] 
C27] 
C28] 
C29] 
C30] 
C31  ] 
C32] 
C33] 


9  KOC  CQ39 
9  KOI.  X 

n  PROGRAM  PERFORMS  K OLMOOORO V-SM1 RNOV 
d  GOODNESS  — OF  — F  X  T  TEST  FOR  LOGNORMAL, 
d  ASSUMPTION  ,  X  MUST  DE  DATA  ARRAY 

jf  ORDER  DATA  AND  FXND  SAMPLE  MEAN  AND  VARIANCE  OF  LOG  DATA 

XS-XC4X] 

MOMGX 

d  PUT  DATA  XN  STANDARD  NORMAL  CASE 
ZM  <GX)-X»)*S 
N«-f  X 

L3JII*.XX  +  1 
-»L4 X  |  (XI>H) 

d  FXND  P  VALUE  FOR  EVERY  DATA  POINTS  AND  PUT  XN  ARRAY 

i-i:x«-n-i 
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